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N00014-78-C-0204  in  support  of  the  U.S.  Naval  Postgraduate  School  research 
project  supported  by  the  Naval  Air  Systems  Command  (Air  370)  and  the  Naval 
Material  Command  (EO/MET).  The  work  was  performed  in  support  of  the  Envi¬ 
ronmental  Physics  Group  at  NPS  under  the  direction  of  Professors  G.  E. 
Schacher  and  K.  L.  Davidson. 
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ABSTRACT 


The  rate  of  production  of  sea  water  droplets  by  bubble  bursting  at 
the  ocean  surface  is  virtually  impossible  to  measure  directly.  Laboratory 
simulations  are  also  difficult.  However,  the  surface  flux  can  be  inferred 


from  measurements  of  the  rate  of  change  of  aerosol  densities  and  changes  in 
the  height  of  the  marine  boundary  mixed  layer.  Data  from  CEWCOM-78  has 
been  analyzed  to  produce  the  aerosol  surface  flux  volume  spectrum  from  0.5 
to  15  pm'radius  at  a  wind  speed  of  8  m/sec.  Using  this  flux  spectrum  and 
equilibrium  aerosol  spectra  from  JASIN,  similar  flux  spectra  are  calculated 
for  sind  speeds  from  5  to  1 5  m/sec. 
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A.  INTRODUCTION 

There  are  two  components  to  aerosols  in  the  marine  boundary  layer: 

1)  continental  (background)  and  2)  locally  generated  sea-spray  droplets. 
The  sea-spray  droplets  are  generated  primarily  by  the  bursting  of  bubbles 
at  the  sea  surface.  The  bubbles  are  produced  by  biological  activity, 
chemical  reaction  and  breaking  waves  (white  caps).  At  wind  speeds  greater 
than  about  3  meters  per  second,  white  caps  are  the  primary  contributor 
to  the  bursting  bubbles.  Thus,  the  sea  surface  is  a  continuous  source  of 
sea-salt  aerosols  in  the  marine  boundary  layer.  These  surface  produced 
aerosols  can  be  characterized  by  a  surface  flux  spectrum,  F  (r),  which 
represents  the  volume  of  aerosol  produced  per  square  centimeter  of  ocean 
surface  each  second  as  a  function  of  aerosol  particle  radius.  This  quantity 
is  a  function  of  wi  nd  speed. 

The  continuous  production  of  sea-salt  aerosols  is  balanced  by  several 
removal  mechanisms.  One  obvious  mechanism  is  the  loss  of  particles  as  they 
fall  back  to  the  surface.  This  settling  under  gravity  is  called  "Stokes 
fallout"  and  is  characterized  by  the  Stokes  velocity.  The  particles  are 
transported  vertically  by  turbulence  in  the  marine  boundary  layer  and  main¬ 
tained  at  a  uniform  density  throughout  the  mixed  layer.  The  growth  of  the 
height  of  the  layer  constitutes  another  loss  mechanism  called  "entrainment." 
The  final  loss  mechanism  is  "rainout,"  which  occurs  when  the  particles 
become  condensation  nuclei  in  the  formation  of  clouds. 

Given  the  surface  flux  spectrum  and  a  parameterization  of  the  removal 
mechanisms,  one  could  predict  evolutions  of  the  aerosol  density  spectrum. 
Unfortunately,  the  surface  flux  spectrum  is  not  known.  However,  we  can 
reverse  the  process--that  is,  use  parameterizations  of  the  removal  processes 
and  evolutions  of  the  aerosol  spectrum--to  obtain  estimates  of  the  surface 
flux  spectrum.  This  .report  describes  the  calculation  of  F  (r)  from  data 
taken  by  the  NPS  Environmental  Physics  Group  during  CEWCOM-78  and  JASIN. 
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B.  AEROSOL  SPECTRUM 


The  first  step  to  analyzing  the  evolution  of  the  aerosol  spectrum  is 
to  remove  those  variations  due  to  changes  in  the  ambient  relative  humidity. 
This  is  done  by  transforming  each  spectrum  to  a  reference  relative  humidity 
(RH  =  80*  or  a  saturation  ratio  S  =  0.8).  We  make  the  usual  assumption 
that  the  particle  radius  a  saturation  S  is  given  by: 

r  =  rQ  g(s)  (1  ) 

where  rQ  is  the  radius  at  S  =  0.8  and  g(s)  =  0.81  exp  [ 0 . 066S / ( 1 . 058  -  S)l. 

If  we  have  a  measured  aerosol  number  density  spectrum,  n(r),  then  the 
transformed  spectrum  at  standard  saturation,  n'(rQ),  is: 

~^r-  =  n‘(rQ)  =  n  (r  g(s))  g(s) 
o 

Since  we  prefer  to  work  with  the  volume  spectrum  V(r) 
corresponding  relationship  is: 

-  V(r  )  =  V(rog(s))/g2(S)  (3) 

We  now  must  separate  the  measured  spectrum  into  its  two  basic  components 

1)  the  background  continental  aerosol  of  nonlocal  origin,  Vc,  which 
is  present  above  and  within  the  marine  layer;  and 

2)  the  sea-salt  aerosol  locally  generated  at  the  sea  surface,  V  , 
which  is  well-mixed  throughout  the  marine  layer  (Z  <  h). 

Thus,  we  assume: 


V(r)  =  Vc(r)  +  Vs(r) 

Z  <  h 

(4a) 

V(r)  =  Vc(r) 

Z  >  h 

(4b) 

(2) 


r3  n(r) ,  the 
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The  continental  aerosol  component  is  represented  with  a  Junge  type 
distribution, 

Vc(r)  =  A/r  (5) 

where  A  is  the  continental  coefficient.  Since  there  is  very  little  sea- 
salt  volume  production  for  r  <  0.3  pm,  we  can  use  the  small  size  range  of 
the  spectrum  to  calculate  the  continental  coefficient  (A1  =  rQ V ' ( r  ) ) . 

This  is  documented  in  Figure  1  where  ensemble  averages  of  volume  spectra 
indicated  that  V 1 ( r  )  is  essentially  independent  of  wind  speed  at  rQ=  0.1  pm. 
At  low  level  wind  speeds,  the  rQ  =  0.1  pm  and  rQ  =  0.3  pm  ensemble  average 
spectra  are  described  quite  accurately  by  Equation  5.  The  claim  that  A1 
is  an  accurate  index  of  continental  influence  is  nicely  validated  by  com¬ 
parison  with  atmospheric  Radon  activity  from  CEVJC0M-78  (Figure  2).  Thus, 
locally  generated  sea-salt  component  can  be  calculated  using: 

V^(rQ)  =  V'(rQ)  -  A7rQ  (6) 

The  remainder  of  this  report  will  deal  exclusively  with  the  volume 
spectra  transformed  to  standard  humidity,  V^(rQ).  However,  to  simplify 
the  notation  we  will  drop  the  prime  and  sub-zero  with  the  understanding  the 
Vs  now  refers  to  the  spectrum  at  RH  =  8073 ,  with  radius,  r,  an  implied 
variabl e. 

The  evolution  of  the  aerosol  spectrum  is  described  by  the  following 
equation : 

”  — |  =  F$  -  (We  +  Ws)  Vs  (7) 
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RH=85% 

JASIN 

20  min  u  ave 

Radius  (^m) 


Ensemble  average  total  aerosol  volume  spectra  from  JASIN. 
The  number  to  the  right  of  the  spectrum  is  the  wind  speed 
category  in  meter/second . 
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where  h  is  the  height  of  the  boundary  layer,  W^the  entrainment  velocity 
and  W  is  the  Stokes  velocity: 

Ws  =  1.57  x  10"2  (1  +  1.2  (|f|y)3)r2  g2(s)  (3) 

where  r  is  in  pm  and  Ws  is  cm/sec.  In  this  equation  we  have  left  out  the 
cloud  formation  removal  mechanism  because  it  acts  on  a  much  longer  time 
seal e. 


C.  FLUX  SPECTRUM  ANALYSIS 


The  CEUCOM- 73  data  have  been  used  to  evaluate  the  surface  flux  term 
of  Equation  7.  Since  the  other  terms  of  the  expression  were  measured,  one 
simply  calculates  Fs  using  a  rearrangement  of  the  terms. 

A  period  was  chosen  from  CEWCOM-78  where  all  data  were  available, 
the  synoptic  conditions  were  fairly  stable,  the  wind  speed  was  reasonably 
constant  and  a  good  mixed  layer  was  present  (Table  I).  There  was  a  slowly 
decreasing  continental  influence  indicative  of  steady  NW  winds  (A1  in  Table 
I  and  Figure  2).  Subsidence  was  assumed  to  be  negligible  so  that  Wg  =  dh/dt 
Thus,  the  final  expression  used  is: 


F 


s 


(~  -  U  )  V 
'dt  s  s 


(9) 


The  scatter  and  uncertainty  of  the  aerosol  and  mixed  layer  depth  mea¬ 
surements  introduces  certain  trade-offs  between  averaging  times  and  statisti 
cal  validity.  For  short  averaging  times,  the  random  scatter  of  the 
individual  terms  is  considerable  leading  to  large  variations  in  the  time 
derivative  terms.  By  increasing  the  averaging  time  to  four  hours,  the  F£ 

"si ynal -to-noise"  becomes  more  reasonable.  A  time  series  of  the  individual 
terms  of  Equation  9  is  shown  in  Figure  3  for  two  different  particle  sizes. 
The  average  surface  production  flux  spectrum  is  shown  in  Figure  4.  This 
flux  spectrum  applies  to  che  average  surface  conditions  for  the  entire  20 
hour  period  (average  wind  speed  was  about  8  m/sec). 
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TABLE  I.  CEWCOM-78  meteorological  and  aerosol  data  for  the 
time  period  analyzed  for  this  report.  The  aerosol 
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radius  particles.  The  dashed  line  is  the  average  contribution 
of  the  surface  production  term  and  should  be  the  sum  of  the  other 
three  terms. 
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Tlic  relative  contributions  of  entrainment,  graviational  fallout  and 
surface  production  arc-  nicely  illustrated  by  defining  an  equivalent  surface 
production  vertical  velocity,  W,,  -  F^/Vr  ,  so  that  in  equilibtium  (dVs/dt  =  0) 

Wf  -  We  +  Ws  (10) 

For  the  average  conditions  found  in  the  analysis  period,  the  entrain¬ 
ment  and  Stokes  terms  were  roughly  equal  at  a  particle  radius  of  5  ,.m 
(Figure  5). 

a.  fll'x  V'ECVbM  unc  sp fed 

di.^a .. o!  the  reasonably  constant  wind  speeds  during  the  20  hour 
C-.tC.yM- 7, v  pc-riuJ,  we  were  able  to  improve  the  statistical  certainty  for  the 

F  calculation  Xj  combining  all  the  data  from  the  period,  assuming  that  the 

flux  woulu  be  a  reasonable  representation  for  U  =  8  m/sec  wind  speed.  Thus, 
we  now  have  the  surface  flux  spectrum  at  a  single  wind  speed.  In  order  to 
esti;  ite  the  flux  at  other  wind  speeds,  we  note  that  the  right  hand  side  of 
Equation  10  is  nearly  independent  of  wind  speed  for  equilibrium  conditions. 
Therefore,  the  flux  at  one  wind  speed  can  be  re1  i  ted  to  the  flux  at  other 
wind  speeds  it  the  equilibrium  volume  spectra  are  known. 

V  (U.) 

MV  ■  FS<U2>  VT-rop-  (11) 

We  have  available  from  JASIN  a  large  set  of  ensemble  averages  of 
aerosol  volume  spectra  at  different  wind  speeds  (Figure  6  and  Table  II). 

It  is  a  simple  matter  to  apply  this  data  to  Equation  11,  using  the 

CEWC0M-73  aerosol  flux  and  equilibrium  spectrum  to  generate  the  surface 
volume  flux  spectra  as  a  function  of  wind  speed  (Figure  7  and  Table  III). 


10 


Velocity  (cm /sec 


THE  BDM  CORPORATION 


Figure  5.  Ensemble  average  contribution  of  entrainment  (W  ),  Stokes  fallout 
(Ws)  and  surface  flux  WF  =  F  /V  assuming  a  stlte  of  dynamic 
equilibrium.  ^  b 
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TABLE  II.  Equilibrium  sea-salt  aerosol -spectra 

in  pm2/cm3  as  a  function  of  radius  (pm) 
and  wind  speed  (m/sec). 
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Figure  7.  Ensemble  average  surface  flux  spectra  deduced  from  Equation  11  and 
Table  II  .  The  number  to  the  right  of  the  spectrum  is  the  wind 
speed  in  meters/second. 
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TABLE  III.  Surface  sea-salt  aerosol  flux, 

F  (r),  (pm2 /cur-/ sec)  as  a  function  of 
particle  radius  (pm)  and  wind  speed 
(m/sec) . 
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E.  NON-EQUILIBRIUM 

For  time  periods  of  a  few  hours,  the  aerosol  spectrum  may  not  be  in  a 
state  of  dynamic  equilibrium.  If  we  rewrite  Equation  7  in  the  following 
form : 


dVs  _  _Fs 

dt  Tp  ’  h 


(12) 


where  the  time  constant, 


=  h/(We  +  Ws) 


(13) 


then  we  see  the  analogy  of  aerosols  and  a  capacitor  charged  by  an  applied 
"voltage,"  Fs,  through  a  "resistance,"  h.  In  this  analogy,  the  "capacitance," 
is  (We  +  Ws)"  1 . 

The  response  time  of  the  aerosol  density  is  a  strong  function  of 
particle  radius  because  Ws  ~  r2.  Values  of  t  for  S  =  0.8  and  h  =  400  m 

r 

are  given  in  Table  IV. 


TABLE  IV.  The  equilibrium  time  constant  for  aerosols 
at  S  =  0.3,  h  =  400  m  and  We  =  0.4  cm/sec. 
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The  boundary  layer  mixing  time,  t  ,  which  represents  the  time  required 
for  changes  in  aerosol  density  to  be  evenly  distributed  throughout  the 
marine  layer  is: 

xm  =  h/W*  (14) 
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"3H 


For  the  analysis  period  (h  =  400  m  and  W*  =  0.6  m/sec)  we  find  =  0.16 
hours.  Therefore,  short  term  variations  in  mixing  volume  (dh/dt)  will 
lead  to  ciianges  in  the  aerosol  density  because  the  production  response  time 
is  much  slower  than  the  mixing  time.  Thus,  for  time  periods  on  the  order 
of  one  hour,  changes  in  the  aerosol  density  (h  dVs/dt)  will  be  highly  cor¬ 
related  with  the  mixing  volume  term  (Vs  dh/dt). 

This  effect,  which  is  particularly  noticeable  for  smaller  particles 
(see  Figure  3  and  Table  III)  is  shown  in  Figure  8.  ^ 

1 
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Figure  8.  Changes  in  aerosol  spectral  density  (dV  / d t )  versus  changes  in 
mixing  volume  (-dh/dts).  This  graph  illustrates  the  dominance 
of  these  terms  in  Equation  9  for  short  time  periods. 
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